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Abstract 

We consider the light scattering from ultracold atoms trapped in an optical lattice inside a cavity. 
In such a system, both the light and atomic motion should be treated in a fully quantum mechanical 
way. The unitary evolution of the light-matter quantum state is shown to demonstrate the non- 
trivial phase dependence, quadratic in the atom number. This is essentially due to the dynamical 
self-consistent nature of the light modes assumed in our model. The collapse of the quantum state 
during the photocounting process is analyzed as well. It corresponds to the measurement-induced 
atom number squeezing. We show that, at the final stage of the state collapse, the shrinking of the 
width of the atom number distribution behaves exponentially in time. This is much faster than the 
square root time dependence, obtained for the initial stage of the state collapse. The exponentially 
fast squeezing appears due to the discrete nature of the atom number distribution. 
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I. INTRODUCTION 



The interaction of quantized light with ultracold quantum gases represents a paradig- 
matic example of the ultimate quantum level of the light matter-interaction, where the 
quantizations of both light and matter (i.e. atomic motion) play equally important roles. In 
particular, a fully quantum description should be applied to describe both the state evolu- 
tion and measurement process. The measurement of photons scattered from ultracold atoms 
provides an example, where the concept of quantum state collapse can be used to describe 
the measurement. Due to the recent experimental developments jl-5], the realization of 
such schemes should become practical. 

In our previous works, we considered the light scattering from atoms trapped in an 
optical lattice into a cavity. We have shown that the detection of photons leaking from the 
cavity represents a quantum nondemolition (QND) measurement of various atomic variables 
observing the light. Such a quantum measurement scheme can be used not only for the 
detection purposes but also for the active quantum state preparation, due to the 
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light-matter entanglement and measurement back-action 

In this paper, we present the details of the atomic state collapse during the photodetection 
process by analyzing the narrowing of the atom number distribution in time. We also 
comment on the short-time unitary evolution of the atomic state before the collapse process 
becomes important. 



II. THEORETICAL MODEL 



We consider the model presented in Refs. [1 lMl3j| : N ultracold atoms in an optical 
lattice of M sites formed by strong off-resonant laser beams. A region of K < M sites is 
illuminated by a weak external probe, which is scattered into a cavity. Alternatively, this 
region is illuminated by the cavity field appearing due to the presence of the probe through 
the cavity mirror. 

We use the open system approach for counting photons leaking the cavity of decay rate 
k. When a photon is detected, the jump operator (the cavity photon annihilation operator 
ax) is applied to the quantum state: \^ c (t)) — > ai\^ c {t)). Between the counts, the system 
evolves with a non-Hermitian Hamiltonian. Such an evolution gives a quantum trajectory 
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for \^ c (t)) conditioned on the detection of photons. 

The expression for the initial motional state of atoms reads 

|*(0)) = ^c g | gi ,.., ?M ), (1) 
i 

which is a superposition of Fock states reflecting all possible classical configurations q = 
{(ft, •-, <Zm} °f N atoms at M sites, where qj is the atom number at the site j. As we 
have shown in Refs. [ill . I12I ]. the solution for conditional wave function takes physically 
transparent form, if the following approximations are used: atomic tunneling is much slower 



than light dynamics and the probe waves are in the coherent state. The tunneling dynamics 
is indeed a very important process in many systems, e.g., in the double- well lattices jl^j ]. 
However, assuming here it is to be negligibly small, one can focus on a different type of the 
system dynamics: the state collapse, which exists even if the tunneling is neglected. The 
conditional state after the time t and m photocounts at t\,t2, ..t m is given by the quantum 
superposition of solutions corresponding to the atomic Fock states in Eq. ([I]): 



\y c (m,t)) = -j^ s ^a q {t l )a q {t 2 )...a q {t 
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where 



xe^c° q \ qi ,...,q M )\a q (t)), (2) 
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dt. (4) 
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If the first photon is detected after the time 1/k, this solution simplifies further |12|: 



*c(rM)) = -^y^<e^W C ;| gi ,...,g M )|a g ), (5) 
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%(t) = ~\a q \ 2 Kt + (rja* - iU 10 a D q 10 a* - c.c.)t/2, (7) 



3 



where a q is the cavity light amplitude corresponding to the classical configuration q (starting 
from here we use the slowly varying amplitudes). It is simply given by the Lorentz function 
(jnj) well-known from classical optics, where a is the external probe amplitude, rj is the 
amplitude of the probe through a mirror; U[ m = gig m /A a (l,m = 0, 1), where gi >0 are the 
atom-light coupling constants, A a = cu\ — u a is the cavity-atom detuning; A p = u p — U\ 
is the probe-cavity detuning. Df m = J2f=i u K T j) u m( r j)Qj are the probe-cavity coupling 
coefficient and dispersive frequency shift that sums contributions from all illuminated atoms 
with prefactors given by the light mode functions tto,i(r). Except the prefactors associated 
with m photodetections a™, the components of quantum superposition in Eq. fl5j) acquires 
the phases contained in d7j). F(t) is the normalization coefficient. 

For several particular cases, the solution for the time- dependent probability distribution 
of atoms corresponding to the state (jSJ) can be simplified further ll|, |l2j. If the probe, 



cavity, and lattice satisfy the condition of the diffraction maximum for light scattering, the 
probability to find the atom number < z < N in the lattice region of K sites is given by 

p(z, m, t) = z 2m e- z2T p (z)/F 2 (m, r), (8) 

z 

with t = 2\C\ 2 K,t, C = iUiodo/(iAp — k), Po(z) is the initial distribution, and F provides the 
normalization. The light amplitude corresponding to the atom number z is a z = Cz. 

If the condition of a diffraction minimum is satisfied, the probability to find the atom 
number difference between the odd and even sites in the lattice region of K — M sites is given 
by the same Eq. (151) . but with a different meaning of the statistical variable —N < z < N . 



III. STATE COLLAPSE FOR SMALL ATOM NUMBER 



Equation (|Sj) shows that any initially broad atom-number distribution po(z) shrinks and 
approaches some very narrow distribution, which corresponds to the atom number squeezing. 
In Refs. 
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121 ]. we have shown that the cental value of the final distribution z$ is given by 



'm/r. Moreover, the full width at half maximum (FWHM) of the distribution can be 
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estimated as 5z w ^/21n2/r, which shows the shrinking of the distribution in time 11 
This type of shrinking (as y/TJr) is rather typical for such photocounting schemes 15]. This 
approximate formula works well if one assumes that (i) the distribution is already rather 



narrow, 8z <C zq, but still (ii) the atom-number probabilities in Eq. flSJ) can be replaced by 
the continuous functions of z, which is a good approximation for very large atom numbers. 
For the large atom numbers, in Refs. jlOy we have demonstrated that the state collapse can 
be very easily described analytically. 

However, it is clear that while the distribution function continues to shrink, its width 
can reach the values of 5z ~ 1, where the approximation of continuous functions obviously 
fails. In this case, one should explicitly take into account the discrete nature of the atomic 
ensemble. If one starts with a macroscopic quantum gas as in Refs. (lOj, it is probably 
not very practical to expect that the distribution p(z, m, t) can really shrink to the widths 
of Sz ~ 1, because various distractive mechanisms will prevent the final state collapse 
to the many-body Fock state \zq,N — zq) with the precisely known atom number z a ^ 
the K illuminated lattice sites and the rest N — zq atoms at the rest of M — K sites. 
Instead, a distribution with some rather small Sz will establish. Thus the approximation 



5z « ^2 ln2/r can work well for the large atom numbers even till the last stage of the 
conditional time evolution. However, if the atom number and number of illuminated lattice 
sites are s,„a„, the situation 1 - reasonable and praetica, elementally Q. 

We have carried out numerical simulations using the Quantum Monte Carlo Wave Func- 



tion (QMCWF) simulations method 
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13) . The numerical results at a single quantum 



trajectory are presented in Fig. 1. They clearly show that the time evolution of the distri- 
bution function width has two different parts. First, the approximation of the square root 
decrease works very well. After that, however, the decrease and, thus the state collapse, 
becomes much faster. In the logarithmic scale figure, it is clear that the shrinking of the 
distribution becomes even exponential. 

Such a behavior of the distribution function can be explained using the fact that the 
discrete functions of z should be taken into account. Let us assume that the atom number z 
changes discretely with a step Z. For example, for the atom number at K sites, which can 
be measured at the directions of diffraction maxima, the step is Z — 1 atom. In contrast, 
for the atom number difference between odd and even sites, which can be measured at the 
directions of diffraction minima, the step is Z = 2 atoms, which reflects the total atom 
number conservation jf], [ill [3]. Note that, when the width of the distribution function is of 
the order of 1, introducing the FWHM becomes meaningless [e.g., the values of p(zo±Z) can 
be already less than p(zo)/2]. Thus, the distribution width should be characterized directly 
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by the square root of the atom number variance: 



(9) 



z 



The insight in the exponential shrinking can be made as follows. If one assumes that at the 
final stage of the state collapse only the atom numbers of zq and zq ± Z are non- negligible, 
the Eq. ([9]) reduces to 



Using Eq. ([8]) (expanding the expressions in Taylor series in Z/zq, taking into account that 
the normalization factor F is a function of m, r, and zq, and substituting m as m = Zqt), one 
can get an estimation (z 2 ) — (z) 2 ~ exp (— Z 2 r). This expression supports the exponential 
shrinking of the atom number distribution width at the final stage of the quantum state 
collapse, which is demonstrated in Fig. 1. 

IV. UNITARY DYNAMICS OF THE QUANTUM STATE 

In this section we address a question about the unitary (coherent) evolution of the light- 
matter state in the short-time limit, where the photon escapes are not important (t < 1/k). 
We will be interested in the scattering into the direction of the diffraction maximum, where 
Df Q = Nx is the fluctuating atom number at K lattice sites. We will consider only the 
transverse probe a (while probing through the mirror does not present, rj = 0). The 
dispersion frequency shift UnD^ will be neglected in this configuration. Using the general 
solution (T51|4]) with those assamtions, and setting k = 0, m = 0, one gets 



(z 2 ) - (z) 2 « Z 2 p(z -Z) + Z 2 p(z + Z). 



(10) 




(11) 



where 



a q (t) = CD q 10 (1 - e lApt ) e- iuJpt 



(12) 



iA p \CD q 10 \H + i\CD 
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sin Apt, 



(13) 
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and C = Uioao/A p . For times such that A p t 3> 1, the phase gets even a simpler form: 
$>q(t) = — iAp\CDf \ 2 t. Thus, the evolution of the state for the interaction at the diffraction 
maximum reads: 



q 

An important property of this solution is that the state components with various atom 
numbers at K lattice sites have different phase evolutions, which quadratically depends on 
the atom number at K sites Nk- For the superfluid state with the large atom number, if 
the number of illuminated sites is much smaller than the total site number, K -C M, the 
stages o<N K is nearh Po.sso.han fl^. Thus, the situation rese.nb.es the p r ohle m s 
of a macroscopic BEC, where the collapse and revival of the matter field can be observed 



171 ] . In both cases, the nontrivial dynamics is a consequence of the particular (i.e. quadratic 
in the atom number) modification of the phase evolution of the atom number components, 
initially constituting the Poissonian atom number distribution. The physical reasons of the 
phase modifications in two cases are however completely different. In the case of Ref. it 
is the atom-atom interaction that brings the nonlinearity in the problem in the form of the 
phase prefactor exp[— iUN(N — l)t/2h], where U is the atom-atom interaction energy. In 
our case, the atom-atom interaction has been completely neglected, and it is the atom-light 
interaction which leads to the phase dependence quadratic in Nk- As can be trace from 
the full solution (@|, the phase depends on the product of probe- mode coupling coefficient 
proportional to D\ Q (which is Nk in this particular case), and the cavity mode amplitude a q . 
As we consider the light amplitudes as dynamical variables (and not as prescribed quantities 
as it is made in many problems of atoms in strong optical lattices), the light amplitudes 
are proportional to the atom number Nk- Thus, the product of the coupling coefficient 
and light amplitude gives us the term quadratic in the atom number Nk- In other words, 
one can say that the dynamical nature of the light mode leads to the "effective" atom-atom 
interaction. 

It is probable, that similarly to Ref. 



171 ]. one can obtain the collapse and revival of the 
matter field in our system, if some matter wave interference measurement will be carried 
out. Since we use the basis with the fixed total atom number N (in contrast to Ref. 
where the coherent atomic state was assumed for the BEC), the calculation of the matter 



interference will be non-trivial. To calculate the interference pattern, one needs to include 
the atom counting procedure. This procedure is well-developed (see, e.g., Refs. flQ) and 
leads to the results rather similar to what is expected from the coherent-state approximation 
for the initial atomic state. Besides that, the methods to disentangle light and matter can 
be used to observe the matter wave interference. 



V. CONCLUSIONS 



We considered the light scattering from ultracold atoms trapped in optical lattice inside 
a cavity. In such a system, both the light and atomic motion should be treated in a fully 
quantum mechanical way. The unitary evolution of the light-matter quantum state was 
shown to demonstrate the non-trivial phase dependence, quadratic in the atom number. 
This is essentially due to the dynamical self-consistent nature of the light modes assumed 
in our model. The collapse of the quantum state during the photocounting process has 
been analyzed as well. We have shown that, at the final stage of the state collapse, the 
shrinking of the width of the atom number distribution behaves exponentially in time. This 
is much faster than the square root time dependence, obtained for the initial stage of the 
state collapse. The exponentially fast atom number squeezing appears due to the discrete 
nature of the atom number distribution. 
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Time r Time r 

FIG. 1: 



Fig.l. Width of the atom number distribution function during the photodetection. De- 
creasing width corresponds to the atom number squeezing. At the initial stage, the shrinking 
is as 1/ y/r, while at the final stage the shrinking is exponential, (a) Quantum trajectory 
without quantum jumps; (b) quantum trajectory with quantum jumps. Photodetection is 
at the diffraction minimum. Total number of atoms N = 100, K = M = 100 sites. 
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